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We present a multiparameter generalization of the Stackel transform (the latter is also known as 
the coupling-constant metamorphosis) and show that under certain conditions this generalized Stackel 
transform preserves Liouville integrability, noncommutative integrability and superintegrability. The 
corresponding transformation for the equations of motion proves to be nothing but a reciprocal trans- 
formation of a special form, and we investigate the properties of this reciprocal transformation. 

Finally, we show that the Hamiltonians of the systems possessing separation curves of apparently 
very different form can be related through a suitably chosen generalized Stackel transform. 
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Introduction 



I 

5^ , 136] for more recent developments, is a powerful tool for producing new Liouville integrable systems from the 
known ones. This is essentially a transformation that sends an n-tupie of functions in involution on a 2n- 
dimensional symplectic manifold into another n-tuple of functions on the same manifold, and these n new 
functions are again in involution. In its original form the Stackel transform affects just one coupling constant 
which enters the Hamiltonian linearly and interchanges this constant with the energy eigenvalue, see [lTlfTB] . 

In the present paper we introduce a multiparameter generalization of the classical Stackel transform, 
which, just like its known counterpart, enables us to generate new Liouville integrable systems from the 
known ones or bring known integrable systems into a simpler form. Unlike the original Stackel transform [TTl 
[18] this multiparameter generalized Stackel transform allows for the Hamiltonians being nonlinear functions 
of several parameters. These properties considerably increase the power of the transform in question. 

Most importantly, under certain natural assumptions the multiparameter generalized Stackel trans- 
form preserves Liouville integrability, superintegrability and noncommutative integrability, see Proposi- 
tions [U and [2] and the discussion thereafter. 
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Moreover, in Section H] we show that the transformations for equations of motion induced by the mul- 
tiparameter generalized Stackel transform are nothing but reciprocal transformations. This generalizes to 
the multiparameter case the earlier results of Hietarinta et al. [H] on the one-parameter Stackel transform. 

The significance of reciprocal transformations in the theory of integrable nonlinear partial differential 
equations is well recognized. These transformations were intensively used in the theory of dispersionless 
(hydrodynamic-type) systems as well as in the theory of soliton systems, see e.g. [291 EI] an d references 
therein. On the other hand, some particular examples of transformations of this kind for finite-dimensional 
Hamiltonian systems are also known, for instance the Jacobi transformation, see [21] and a recent survey 
|36j . The reciprocal transformations of somewhat different kind have also appeared in [T8| 138 } [35]. 

In the present paper we consider reciprocal transformations for the Liouville integrable Hamiltonian 
systems in conjunction with the generalized Stackel transform and, in contrast with the earlier work on the 
subject, we concentrate on the multi-time version of these transformations. 

In fact, as we show in Section H] below, these transformations, when applied to the equations of motion 
of the source system, in general do not yield the equations of motion for the target system unless we 
restrict the equations of motion onto the common level surface of the corresponding Hamiltonians, see 
Propositions [3] and H] below for details. 

We further show that for two Liouville integrable systems related by an appropriate multiparameter 
generalized Stackel transform for the constants of motion we have the reciprocal transformation relating 
the corresponding equations of motion restricted to appropriate Lagrangian submanifolds, see e.g. Ch.3 of 
[12J and references therein for more details on the latter. 

Moreover, we present a multitime extension of the original reciprocal transformation from [T5] . and 
study the applications of this extended transformation to the integration of equations of motion in the 
Hamilton- Jacobi formalism using the separation of variables, cf. [TT] . 

In the rest of the paper we consider the relations among classical Liouville integrable Stackel systems on 
2n-dimensional phase space. In [7J infinitely many classes of the Stackel systems related to the so-called seed 
class, namely, the fc-hole deformations of the latter, were constructed. Here we show that any /c-hole defor- 
mation can be obtained from the Benenti-type system through a suitably chosen multiparameter generalized 
Stackel transform, and present the explicit form of the transform in question along with its inverse. 

2 Multiparameter generalized Stackel transform: 
definition and duality 

Let (M,P) be a Poisson manifold with the Poisson bracket {/, g} = (df,Pdg). Consider r functionally 
independent Hamiltonians Hi, i — 1, . . . ,r, on M, and assume that these Hamiltonians further depend on 
k < r parameters «i, . . . , so 



where x G M. Note that in general r is not related in any way to the dimension of M except for the obvious 
restriction r < dimM; see, however, the discussion after Proposition [TJ Also, in what follows all functions 
will be tacitly assumed to be smooth (of the C°° class). 
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Suppose that there exists a fc-tuple of pairwise distinct numbers Sj G {1, . . . , r} such that 



det(|3ff^0a^=i,... lfc )^O. (2) 

Now fix a fc-tuple {si, . . . , s k } such that ([2]) holds and consider the system 

H Si (x, aci, . . . ,a k ) = a h z = 1, . . . , k, 

where fij are arbitrary parameters, as a system of algebraic equations for ai, . . . ,a k . By the implicit function 
theorem, the condition (j2J) guarantees that the solution of this system exists and is (locally) unique. We 
can write this solution in the form 

a.i = Ai(x, dii, ■ ■ ■ , a/c), i — l,...,k. 

Now define the new Hamiltonians H Si , i = 1, . . . , k, by setting 

H 8i = Ai(x, 5i!,. . . ,a k ), i = l,...,k. 

In other words, the Hamiltonians H Si , i — 1, . . . , k are defined by means of the relations 

H Si \[<s>) = on, i = l,...,k. (3) 

Here and below the subscript [<&] means that we have substituted H Sj for oti for all % = 1, . . . , 
Next, let 

Hi = Hi\[$\, i = l,...,r, Sj for j = 1, . . . , fe. (4) 

Note that the Hamiltonians Hj involve k parameters aj, i = 1, . . . , k for all j = 1, . . . , r: 

Fj = Hi(x,a u . . . , fife), 2 = l,...,r. 

We shall refer to the above transformation from ifj, z = 1, . . . ,r, to ifj, z = 1, . . . ,r, as to the k- 
parameter generalized Stackel transform generated by H Sl , . . . , H Sk . In analogy with [11] we shall say that 
the r-tuples Hi, i — 1, . . . , r, and Hi, i = 1, . . . , r, are Stackel- equivalent. 

The condition (T5]) guarantees that the above transformation is invertible. Indeed, consider the dual of 
the identity that is, 

H St \y^=ai, i = l,...,k, (5) 

where the subscript [<3>] means that we have substituted H s . for fij for all i — 1, . . . , k. 

Moreover, the functional independence of the original Hamiltonians Hi, i = 1, . . . ,r, implies the func- 
tional independence of Hi, i = 1, . . . , r. Indeed, the functional independence of H, % = 1, . . . , r, means that 
dim span (dffj, i = 1, . . . , r) = r on another open dense subset {7 of M. Using (|2J), © and (j4j) we readily 
see that this implies dimspan(<iifj, i = 1, . . . , r) = r on another open dense subset £/ C U of M. In turn, 
the latter equality means nothing but the functional independence of Hi, % = 1, . . . , r we sought for. 

Let us stress that here and below the differentials are computed under the assumption that the parame- 
ters are considered to be constant, i.e., if H = H(x, ot\, . . . , a k ) then in the local coordinates x b on M we have 

dimM f)H 
6=1 
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By the implicit function theorem the condition (j2J) guarantees that we can solve (jSJ) with respect to H Sj , 
j — 1, . . . , k. If we do this and define the remaining Hamiltonians Hi by the formulas 



Ha = Hj 



[*]' 



i = 1, • • • ,r, i / for j = 1, k, 



(6) 



then it is straightforward to verify that and (jl]) hold identically In other words, the formulas (jSJ) and 
(J6]) define the inverse of the transformation defined using ([3]) and (HI). 

Clearly, these two transformations are dual, with the duality transformation swapping Hi and Hi for all 
i = 1, . . . , r and swapping ctj and a, for all j = 1, . . . , fc. 

Note that in the special case when the Hamiltonians Hi are linear in the parameters acj, the above 
formulas undergo considerable simplification, and we can explicitly express Hi via Hi. 

Namely, let 



-i , z = l,...,r. 



Then equations ([3]) take the form 



(7) 



«i, z = 1 



)•••)"') 



and we can readily solve them for i7 s .: 



F s . = detWi/detW, 



(9) 



where W is a k x matrix of the form 



h-W ... H {k) 



H 



(i) 



(fc) 



and PVj are obtained from W by replacing H^ by iij^ — aj for all j = 1, . . . , fc. 
By (jl]) we have 



3=1 



i i 



i = l,...,r, i^s, for j = 1, . . . , fc, 



(10) 



where are given by (Q. It is straightforward to verify that if we set = 1 then the transformation 
given by (Q and ffTUl) becomes nothing but the standard Stackel transform pj] , also known as the coupling- 
constant metamorphosis [IB] . 



3 Multiparameter generalized Stackel transform 
and (super) integrability 

It turns out that the /c-parametric generalized Stackel transform preserves the commutativity of the Hamil- 
tonians Hi. More precisely, we have the following result: 
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Proposition 1 Let H, L , i = l,...,r, be functionally independent and let H i; i = l,...,r, be related to 
Hi, i = 1, . . . , r, by a k -parameter generalized Stdckel transform |3j) ; generated by H S1 , . . . , H Sk , where 
k < corankP + (1/2) rank P. 

Then the following assertions hold: 

i) if {H Si , H Sj } = for all i, j = 1, . . . , k then {H Si , H Sj } = for all i,j = 1, . . . , k; 

ii) suppose that i) holds and for a j G {l,...,r} ; j ^ Si,...,Sk, we have {H s .,Hj } = for all 
i = 1, . . . , k; then {H Si , Hj } = for all i = 1, . . . , k; 

Hi) suppose that for a natural m < corankP + (l/2)rankP we have an m-tuple of pairwise distinct 
integers h, . . . ,l m G {1, . . . , r} such that s q G {h, . . . , l m } for all q = 1, . . . , k, and {H^, H^} = for 
all i, j = 1, . . . , m; then {H^, H^} = for all i, j = 1, . . . , m. 

Before we proceed with the proof of Proposition [TJ some remarks are in order. First of all, corank P + 
(1/2) rankP is easily seen to be the maximal possible number of functions in involution on M with respect 
to the Poisson bracket associated with P. 

Next, from Proposition [1] it is immediate that the transformation defined by ([3]) and (j4]) preserves 
(super) integrability. Namely, under the assumptions of Proposition [U iii) let dimM = 2n, rankP = 2n, 
and m = n. Then the dynamical system associated with any of H^ is Liouville integrable, as it has n 
commuting functionally independent integrals, H^, j = l,...,n, in involution. By Proposition [TJ iii) the 
dynamical system associated with any of Hi. enjoys the same property, the required integrals of motion in 
involution now being Hi t , i = 1, . . . , n. 

Note that if under the assumptions of Proposition [U iii) we have dim M = 2n, rank P = 2n, and 
m < n then, under some technical assumptions and in a suitable vicinity U C M, for the dynamical system 
associated with any of Hi t , i = 1, . . . ,m, there exists a symplectic submanifold fibred into m-dimensional 
invariant tori [27J [281 E] • The tori in question are intersections of this symplectic submanifold with the 
common level surfaces of Hi v i = 1, . . . ,m. Proposition [TJ iii) implies that this property is preserved by 
the multiparameter generalized Stackel transform defined by ([3]) and (HI), i.e., for the dynamical system 
associated with any of H^, i — 1, . . . ,m, there exists, again under certain technical assumptions and in a 
suitable vicinity U C M, a symplectic submanifold fibred into m-dimensional invariant tori. 

Now let dimM = 2n, rankP = 2n, r > n, and suppose that {H Si , Hj} = for all i = 1, . . . , k and 
for all j = 1, . . . , r. Then the Hamiltonian H s . is superintegrable for any j G {1, . . . , k} as it has r > n 
integrals of motion Hi, i = 1, . . . ,r, and by Proposition [lj ii) the Hamiltonian H s . is superintegrable for 
any j G {1, . . . , k} as well, the integrals of motion now being Hi, % — 1, . . . , r. 

Moreover, the multiparameter generalized Stackel transform defined by ([3]) and (jlj) also preserves non- 
commutative integrability in the sense of [261 [9]. We start with the following result: 

Proposition 2 Under the assumptions of Proposition^ iii) suppose that dimM = 2n, P is nondegenerate 
(rankP = 2n), and the algebra T of functions on M generated by H 1} . . . , H r is closed under the Poisson 
bracket and is complete in the sense of [9]. Further suppose that ker{, }|^- = Tq, where is the algebra of 
functions on M generated by H^ , . . . , H\ m . 

Then the algebra T of functions on M generated by H 1 , . . . , H r is also closed under the Poisson bracket 
and complete. 
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Consider an algebra A of functions on a symplectic manifold M and assume that A is closed under the 
Poisson bracket. Recall (see [9] for precise definitions and further details) that the differential dimension 
ddim A of A is, roughly speaking, the number of functionally independent generators of A. The differential 
index dind*4 can be (informally) defined as dind .A = ddimker{, }|_4, and A is said to be complete [H] if 
ddim A + dind A = dim M on an open dense subset U C M. 

Sketch of proof of Proposition^ First of all, it is immediate that the algebra T generated by H\, . . . , H r is 
also closed under the Poisson bracket. As we have already noticed in Section [21 the functional independence 
of Hi, % = 1, . . . , r, implies that of Hi, i = 1, . . . , r, and hence we have ddim T = ddim T = r. In turn, as 
ker{, }|jr = JF , we have dindJF = ddimjF = m. 

By Proposition [H iii) we have {Hi-,Hi.} = 0, i,j = l,...,m, so kerj,}^ D JF , where JF is the 
algebra of functions on M generated by H^, . . . , H\ m . Therefore dindJF > ddimjF = m. However, as 
we obviously have ddim T + dind T < dim M and, on the other hand, we know from the above that 
ddim T + dind T > r + m = dim M, we conclude that ddim T + dind T = dim M, and thus the algebra T 
is indeed complete. □ 

Therefore, if under the assumptions of Proposition [2J there exists an integer i® £ {1, . . . ,r} such that 
{H iQ ,Hj} = 0, j = 1, . . . ,r, and thus the dynamical system associated with H io is completely integrable 
in the noncommutative sense [251 E]? as this system possesses a complete algebra of integrals of motion, 
then so does the dynamical system associated with H io . 

Proof of Proposition^ Prove i) first. For any smooth functions / and g on M that further depend on 
the parameters ai, . . . , we have the following easy identities: 

k 

{/![*]> = if'9}\ m + J2( d f/ da ^w{H Sj ,g}\mi (ii) 

{/![*]> = + E(0//0«j)l[*]{-^->0}l[*] + E ( d 9/daj) I [$]{/, ^ S J|[0] 

(12) 

Using the assumption {H Si , H Sj } = and ([3]), we find that 

= {ai — H s . , atj — H Sj } = {H s . \ [$] — H s . , H S:j \ [$] — } , 

whence 

I [*] - H Si , H Sj | [ ] - H Sj }| [$] =0. 

Writing out the Poisson bracket on the left-hand side of the latter identity using ( flTl) for the brackets 
{if s . | [$] , # Sj . } | m and {H a ., H aj \[$\}\[$\ and ([T2j) for the bracket ^|[$]} we obtain 

^ (aiy Si /^)| m (^ s ./^)| [4] {# Sp ,^ s j| m = o, 

p,q=l 

whence using ([2]) we readily find that for all p, q = 1 , . . . , k we have 

{F Sp ,# s J| m = 0. 
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However, H Sp are independent of a q for all q — 1, . . . , k, so 

{H Si , H Sj } = {H Si , H Sj }\ m = 0, 

and the result follows. 

As we have already proved i), to prove ii) we only need to show that if {H Si , Hj } — for alH = 1, . . . , k 
then {H Si , Hj } = for alii — 1, ... ,k. 

As Hi, i — 1, . . . , r, are independent of a p for all p = 1, . . . , k by construction, we have 

{H St , F J0 } = {H Si , if io } | [$]. 

Moreover, as j ^ s p for all p — 1, . . . , k by assumption, by virtue of (jlj) the relation {H s ., if J0 } | [$] = 
is equivalent to 

{ H si,Hj 1[*]}|[*] = o. 

In turn, using ffTTT) we can rewrite the Poisson bracket {if Si , -Hj 1 [*] } I [*] as follows: 

k 

{H Si ,H jo |[$]}|[$] = {if Si , -Hj }|[$] - 2j(9iI io /9Q!p)|[$] {i^ , if Si } 

p=i 

As {if s , H s .} = by i), we see that 

{if Si , H jo | [$] } | [$] = {if Si , if io }![$]. 
Now, in analogy with the proof of i), consider the identity 

= {a p ,H jo }\ m = {H Sp \ m , H jo }\ m . 
Using ( fTTl) and our assumptions yields 

k 

= i H s p \m, H jo }\m = {dH Sp /dai)\ [$ ] { if Si , g jo } | . 
Finally, using ([2]) we conclude that 

{H 8i ,H jo }\ m = 0, (13) 

whence {H s ., Hj = 0, and the result follows. 

Part iii) is proved in analogy with ii). Namely, in view of i) and ii) we only need to prove that the 
conditions {Hi., H tj } — 0, i, j — 1, . . . , m imply {H[., H^} = for all i, j = 1, . . . , m such that Zj ^ s p and 
lj 7^ s p for all p = 1, . . . , k. 

If U 7^ s p and lj ^ s p for all p = 1, . . . , k then we have 

{H h ,H lj } = {H h \ m ,H lj \ m }. 
Using f|T2|) and f|T3|) for j = li and jo = lj we readily find that 

{Hiilm, Hijlm} = °> 

and the result follows. □ 

Note that the computations in the above proof bear considerable resemblance to those in the theory of 
Hamiltonian systems with second-class constraints, see e.g. the classical book of Dirac [H] . 
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4 Reciprocal transformations for the equations of motion 



Recall that the equations of motion associated with a Hamiltonian H and a Poisson structure P on M read 
(see e.g. [3]) 

dx b /dt H = (X H ) b , b = l,...,dimM, (14) 

where x b are local coordinates on M, Xh = PdH is the Hamiltonian vector field associated with H, and 
ts is the corresponding evolution parameter (time). 

Throughout the rest of this section we tacitly assume that Hi, i — 1, . . . , r, are related to i — 1, . . . , r, 
through the fc-parameter Stackel transform ([3]), (jlj) generated by H Sl , . . . , H Sk . 

Suppose that {H Si , H Sj } = for all i, j = 1, . . . , k ,and consider simultaneously the equations of motion 
( [141) for the Hamiltonians H Si with the times t Si and for H Si with the times t Si : 

dx b /dt s . = (X Hs .) b , b = l,...,dimM, i = l,...,k, (15) 

dx b /dt s . = {X & ) b , b = l,...,dimM, i = 1, . . . , k. (16) 

In analogy with [18] consider a reciprocal transformation (see e.g. [29l [3T| [32] for general information on 
such transformations) relating the times t Si and t Sj : 




Proposition 3 Suppose that k < corankP + (1/2) rankP and {H Si , H Sj } = for all i,j = 1, . . . , k, and 
consider the equations of motion < T75]) for H Si , i — 1, . . . , k, restricted onto the common level surface N & of 
H Si , where 

N & = {x e M\H Si (x,a u ...,a k )=ai, i = 1, . . . , k}. 

Then the transformation ( fj?| ) is well defined on these restricted equations of motion and sends them 
into the equations of motion ( TJ7| ) for H Si , i = 1, . . . , k, restricted onto the common level surface N a of H Si , 
where 

N a = {x e M\H ai (x, ax, ...,a k )=a u i = 1, . . . , k}. 



Note that the level surfaces in question, N a and N & , represent the same submanifold of M, i.e., N a = N & . 
This is readily verified using the relations ([3]) and (JSJ). 

Proof. First of all show that (1171) is well-defined, that is, we have 



d 2 L 



d 2 L 



p,q 



, k, 



by virtue of equations (jToj) restricted onto N & . 
Using (fT?]) we find that ffTS]) boils down to 



OOii 

~dt7 



\ 



I Q I 



doti 



k. 



(19) 
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In turn, using (fT5l) we readily find that (fT9l) takes the form 



dH s 

doti 



No 



dH s 
da; 



[*] 



3 i>p 



iVn 



and the latter equality can be proved by taking the partial derivative of the relation {H Sp ,H Sq } 
respect to c^. 

Next, Eq.flHD yields 



d 

dt.. 



E 



3=1 v J 



d 

dt s . 



i — 1, . . . , k. 



Taking into account (1151) and (1161) we conclude that we have to prove that 

k 



where Ijv- denotes restriction onto N & . 



da i 



x h s .\n & , i — 1, • • • , k, 



As Xjj = PdH for any smooth function H on M, Eq. (120|) boils down to 



[*] 



0, z = 1, . . . , k. 



On the other hand, taking the differential of ([3]) we obtain 



(dH Si )\ m +J2 



OH, 



v da,- 

3=1 x J 



(dF s J|[ $] =0, i = l,...,k. 



As iJ s . are independent of a p , for all p — 1, . . . , k we have (dH Sj )\[<t>] = dH s ., so ( 1221) yields 



E 



5=1 J 



dH Sj = -(dH Si )\ m , 



and fl20l) takes the form 



(P(dH Si -(dH Si )\ m ))\ N _ =0, i = l,...,k. 



In the local coordinates x b on M we have 



(p (dif Si -(d# Si )y) 



dim A/ 



6=1 



6=1 



[*L 



(Pdx b ) 



z — 1, . . . , k. 



By virtue of ([3]) and (jSJ) A/q, and represent the same submanifold of M, whence 



dH,. fdH 



dx b V dx b 



dH,. fdH, 



dx b \ dx b 



dH Si 
dx b 



dH s 

dx b 
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We used here an easy identity 



' 9H Si 
dx b 



dx b 



Thus, the left-hand side of (123]) . and therefore that of fl2T]) . vanishes, and the result follows. □ 
Now assume that all if, are in involution: 



{H l ,H 1 }=0, 
Then by Proposition [TJiii) so are Hi, i.e., 

{H t ,H 3 } = 0, 



*,3 



*,3 



r. 



and we can consider two sets of simultaneous evolutions, 



dx b /dU = {X Hi )\ 

dx b /dti = (X^f, 
and the following extension of f[T7|) : 



dimM, 
dimM, 



(24) 
(25) 



dt s 

tn ■■ 



-E 

3=1 



"q, 



dJL, 
don 

1,2,. 



dt 



[*] 



k. 



(26) 



k. 



r, q 7^ s p for any p = 1, 

In analogy with Proposition [3] we can prove the following result. 

Proposition 4 Suppose that {Hi,Hj} = for all i,j = l,...,r and r < corankP + (l/2)rankP ; and 
consider the equations of motion [Efy for Hi, i — 1, . . . , r, restricted onto N & - 

Then the transformation j[2h)) is well defined on these restricted equations of motion and sends them 
into the equations of motion (25\) for Hi, i = 1, . . . , r, restricted onto N a . 

Note that the transformations from Propositions [3] and H] do not change the dynamical variables x. 
In particular, under the assumptions of Proposition [3] for any given i from 1 to k the trajectories of the 
dynamical system associated with H s . are identical to those of the dynamical system associated with H s ., if 
we consider the trajectories as non-parametrized curves. In other words, the transformation ( flTI) amounts to 
the reparametrization of the times associated with H s . for all j = 1, . . . , k. Notice, however, that in general 
the reparametrization in question is different for different trajectories, as one can readily infer from (Tl7|) . 

As a final remark note that it could be interesting to compare the above reparametrization results with 
those arising in the theory of projectively equivalent metrics [TU1 134"]. 



5 Canonical Poisson structure 

In this section we tacitly assume that Hi, i = 1, . . . ,r, are related to Hi, i = 1, . . . ,r, through the k- 
parameter Stackel transform (0J generated by H S1 , . . . , H Sk . We further assume that M = M, 2n , P is 
a canonical Poisson structure on M, and Aj, /i,, i = l,...,n, are the Darboux coordinates for P, i.e., 
{Xi, fij} = 5ij. Let A = (Ai, . . . , A n ) and n = (fi±, . . . , fin). Then the Hamilton-Jacobi equations for Hi and 
Hi have a common solution, cf . [TTJ . Namely, we have the following generalization of the results of [11] to 
the case of multiparameter generalized Stackel transform: 
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Proposition 5 Suppose that {H Si , H s ,} = for all i,j = 1, . . . , k. Let S = S(X, a\, . . . , a k , E Sl , . . . , E Sk , 
di, . . . , a n _ k ), where a« are arbitrary constants, be a complete integral of the stationary Hamilton- J acobi 
equation for the Hamiltonians H Si = H St (X, fi,a\, . . . , ctk), 

H s . (A, dS/dX, ,a k ) = E s ., i = l,...,k. 

If we set E Sj = cti and at = E Si for all i = 1, . . . , k then S also is a complete integral of the stationary 
Hamilton- J acobi equation for the Hamiltonians H s . = H S .(X, fi,&i, . . . ctk), 

H Si (X,dS/dX,a h ...a k ) = E Si . 

Further assume that r <n, and {Hi, Hj} — 0, i, j — 1, . . . ,r, and 

S = S(X,cti, . . . , a k , Ei, . . . , E r 

) (27) 

where ai are arbitrary constants, be a complete integral for the system of stationary Hamilton- J acobi equa- 
tions 

Hi(X, dS/dX, a 1 ,...,a k )=E i , i = 1, . . . , r. 

// we set 

ctj = E Sj , E Sj = dij, j = l,...,k, and E t = E h i = l,...,r, i ^ s p for all p=l,...,k, 



then S [21) is also a complete integral for the system 

Hi(X, dS/dX, oti,...,a k )=Ei, i = 1,. . .,r. 

This result suggests that the multiparametric generalized Stackel transform potentially is a very powerful 
tool for solving the Hamilton- Jacobi equations (and hence the equations of motion) for Hamiltonian dynam- 
ical systems. Indeed, if we can solve the stationary Hamilton-Jacobi equations for the original Hamiltonians 
Hi, then by Proposition [5] we can do this for the transformed Hamiltonians Hi as well, and vice versa. 

As for the equations of motion, in addition to general Propositions [3] and IH a somewhat more explicit 
result can be obtained by straightforward computation: 

Corollary 1 Suppose that r = n, {Hi, Hj} = for all i, j = 1, . . . , n, d 2 Hi/dajdfi = for all i — 1, . . . ,n 
and all j = l,...,k, and that Xj, j = l,...,n, can be chosen as local coordinates on the Lagrangian 
submanifold Ne = {(A, /j,) G M|fZj(A, fi, a%, . . . , a k ) = E t , i = l,...,n} (in other words, the system 
Hi(X, fi, , a k ) = Ei, i = 1, . . . , n, can be solved for [i ), and that we have 

ctj = E Sj , E Sj = ctj, j = 1, . . . , k, and Ei = Ei, i = 1, . . . , n, i ^ s p for all p = 1, . . . , k. (28) 

Then the reciprocal transformation If2h]) turns the system 

dX/dU= (dHi/dfj,)\ NE , i = l,...,n, (29) 

into 



dX/dti = (dHi/dfj,)\^_, i = l,...,n, (30) 
where N E = {(A,^t) G M\Hi{X, fx, ax, . . . ,a k ) = E { , i = 1, . . . ,n}. 
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Recall that Ne and N E in fact represent the same Lagrangian submanifold of M, cf. the remark after 
Proposition [3j 

For instance, if we have k = 1, a% = a, si — s, and take 

H = Gi(X)ii) + Vi(X) + aWi(\), i = l,...,n, (31) 

where (-, •) stands for the standard scalar product in M n and Gi(X) are n x n matrices, then the system 
(|29|) reads 

dA/cft* = Gi(X)M, (32) 

where /x = M(X, a, Ei, . . . , E n ) is a general solution of the system Hi(a, A, n) = E iy % = 1, . . . , n. 
If we eliminate M from fl32|) then we obtain the dispersionless Killing systems (cf. [15l [16]) 

A ti =G i (G s y 1 \ t3 , i = l,2,...,s-l,s + l,...,n, (33) 

and the reciprocal transformation (1261) . which in our case reads 

n 

tit s = - ) j Wi(\)dti, U = U, i ^ s, 
i=i 

turns (1331 into 

A t - = GiiG.)- 1 ^, i = 1, 2, . . . , 8 - 1, 8 + 1, . . . , n, (34) 

where the quantities G s = — G s /W s and Gi = Gi — WiG s /W s , i — 1, 2, . . . , s — 1, s + 1, . . . , n, are related to 
the Hamiltonians 

Hi = ^(/x, Gi(\)n) + Vi(X) +&Wi(\), i = l,...,n, (35) 

which are Stackel-equivalent to H iy i = 1, . . . , n. 

We can now apply Proposition [5] in order to obtain the solutions of equations of motion (129|) and (130]) 
as follows: 

Corollary 2 Under the assumptions of Corollary^ suppose that 

S = S(X, ax,... ,a k ,E x , ... ,E n ) (36) 

is a complete integral for the system of stationary Hamilton- J acobi equations 

H(X, dS/dX, ai, . . . , oik) = E h i = 1, . . . , n. 

Then a general solution of (HP)] for i = d can be written in implicit form as 

dS/dE 3 =5 jd t d + b ]} j = l,...,n, (37) 

where bj are arbitrary constants, and by virtue of a general solution of for i = d can be written 
in implicit form as 

dS/OEj = 5 jd t d + bj , j = 1 , . . . , n. (38) 
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Comparing (j37j) and fl38|) and using fl28l) we readily see that, in perfect agreement with (1261) . = for 
i ^ Si, . . . , s k , but £ s = dS/dE Sj — 6 S = dS/daj — b s . while t Sj = 3S/dE s . — 6 Sj = dS/daj — 6 Sj . Thus, 
the above approach does not yield an explicit formula expressing t Sj as functions of A, /j,, and t Si . 

In order to find a complete integral (1561) we can use separation of variables as follows (see e.g. J33J CTJ an d 
references therein). Under the assumptions of Corollary [2] suppose that Aj, fa, i = 1, . . . ,n, are separation 
coordinates for the Hamiltonians Hi, i = 1, . . . , n, that is, the system of equations Hi(\, (/,, ax, . . . , a k ) = Ei, 
i = 1, . . . , n, is equivalent to the following one: 

ipi(Xi, fa, ai,..., a k , Ex,..., E n ) = 0, i = l,...,n, (39) 

which is nothing but the set of the separation relations^ on the Lagrangian submanifold iV^. 
On the other hand, under the identification (|28|) the system (139|) is equivalent to 

Hi{\, //,ai,..., oc k ) = E h i = l,...,n. (40) 



Thus, the Stackel-equivalent n-tuples of Hamiltonians share the separation relations f!39|) provided fl28|) 
holds. 

Consider the system of stationary Hamilton-Jacobi equations for Hi 

Hi(\,dS/d\, a u . . . ,a k ) = E h i = l,...,n. (41) 
By the above, (j4T|) is equivalent to the system 

<fi(\i,dS/dXi,ai, . . . ,a k ,E h . . . ,E n ) = 0, i = l,...,n. (42) 
Suppose that (!39|) can be solved for fa, i = 1, . . . , n: 

fa = Mi(Xi, ai,..., a k , Ei,..., E n ), i = l,...,n. 
Then there exists a separated complete integral of ( l4*2l) . and hence of (1411) . of the form (cf. e.g. [7]) 



n „ 

S = *y] Mi(X h ai,..., a k , Ei,..., E n )dXi 
i=i J 



(43) 



and general solutions for ( 1291) and ( 13U|) can be found using the method of Corollary [2j 
In this case the formulas (1371) take the form 



n „ 

/ (9Mi(Xi, ai,..., a k , Ei,..., E n )/dEj)dXi = 5 jd t d + b h j = l,...,n, (44) 
i=i J 

and expressing Xi as functions of t& from (T441) is nothing but an instance of the Jacobi inversion problem. 
In particular, for d = we have 



n „ 

f Sj + b Si = dS/dE Si = dS/dai = J2 W(A,, «!,..., a fc , ^, . . . , E^/da^dh 

i=i J 



i = l k. 



^^Note that the separation relations involving parameters appear, in a rather different context, in the paper [37j where they 
are employed for the construction of separation variables. 
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6 Multiparameter generalized Stackel transform 
and deformations of separation curves 



Under the assumptions of Corollary [fl suppose that Aj, fa, i = 1, . . . ,n, are separation coordinates for the 
n-tuple of commuting Hamiltonians Hi, i = 1, . . . , n. Then the Lagrangian submanifold Ne is defined by n 
separation relations (159]) . Further assume that all functions </?j are identical: 

ipi = (p(\ i ,Hi,ai,...,ak,E 1 ,...,E n ), i = 1, ...,n. (45) 

Then relations fl39l) mean that the points (Aj, i — 1, . . . , n, belong to the separation curve [33], [7J 

cp(X, ft, , a fc , £„) = 0. (46) 

If the relations 

tp(\i, fii, ai,..., a k , Hi, ... , H n ) = 0, i = l,...,n, 

uniquely determine the Hamiltonians Hi for i — 1, . . . , n, then for the sake of brevity we shall say that H^ 
for z = 1, . . . , n have the separation curve 

(p(X, fj,,a h ..., a k , H u . . . , H n ) = 0. (47) 

Fixing values of all Hamiltonians Hi = E^, i = l,...,n, picks a particular Lagrangian submanifold from 
the Lagrangian foliation. It is also clear that the Stackel-equivalent n-tuples of the Hamiltonians Hi, 
i = 1, . . . , n, and Hi, i = 1, . . . , n, share the separation curve (14TI) provided ([3]) and ([5]) hold. 

In the rest of this section we shall deal with a special class of separation curves of the form (cf. e.g. [7] 
and references therein) 

n 

Y,H^=^X,fi), (48) 
i=i 

where (3j are arbitrary pairwise distinct non-negative integers, fti > fa > • • • > ft n . In fact one always can 
impose the normalization fa = by dividing the left- and right-hand side of fj48|) by if necessary, but 
we shall not impose this normalization in the present paper. 

For a given n, each class of systems fT48|) is labelled by a sequence (fa, ... , fa) while a particular system 
from a class is given by a particular choice of ip(\,n). In particular, the choice ip(X,jji) = |/(A)/i 2 + 7(A) 
yields the well-known classical Stackel systems. All these systems admit the separation of variables in the 
same coordinates (A,,/ij) by construction. 

We shall refer to the class with the separation curve 

n 

^2 Hj^ nj — V'(^) y) (49) 

as to the seed class. Note that if ^(A, /i) = |/(A)/i 2 + 7(A) we obtain precisely the Benenti class of Stackel 
systems [HE]- The seed class is a rather general one: it includes the majority of known integrable systems 
with natural Hamiltonians [7j. 

It turns out that, roughly speaking, the n-tuple of Hamiltonians having the general separation curve 
(T4H]) can be related via a suitably chosen generalized multiparameter Stackel transform to an n-tuple of 
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Hamiltonians having the separation curve (j49H from the seed class. The exact picture is a bit more involved, 
as in fact we need to consider the deformations of the curves in question. 
Define first an operator R[ that acts as follows: 



R[(F) = F + f\ k - (X k /k\){d k F/dX k 



A=0- 



For instance, we have 



Now let 



j=0J¥=k 



n _ n _ 

F = E Hj\ n ~ j and F = £ HjX n ~ j . 

3=1 3=1 



For any integer m define [7] the so-called basic separable potentials K by means of the relations 



(m) xn~j 







(50) 



that must hold for A = Aj, i = 1, . . . , n. 

Under the assumptions of Corollary [TJ consider an n-tuple of commuting Hamiltonians of the form 

k 



3=1 

where 7^, j — 1, . . . , k, are pairwise distinct integers. 

Suppose that the Hamiltonians (15"T1) have the separation curve of the form 

k 

EM 7j ' + f o = ^(A,/i), 



(51) 



(52) 



where 7, > n — 1 for all j = 1, . . . , k, and ji 7^ jj if i 7^ j for all i, j = 1, . . . , k. 

Now pick k < n distinct numbers Sj G {l,...,n} and define the Hamiltonians Hi by means of the 
following separation curve 



E ^ A^ + Ki Sl • • ■ iC Sfc (F ) = ^(A, n). 
3=1 



(53) 



This means that Hi are the solutions of the system of linear algebraic equations obtained from ( 1531) upon 
substituting A.; for A and Uj for u into (153]) for i = 1, . . . , n. 

Proposition 6 Under the above assumptions the n-tuple of Hamiltonians Hi, i = l,...,n, is Stackel- 
equivalent to H i; i = 1, . . . , n. 



The n-parameter generalized Stackel transform relating Hi, i 



n to Hi, i — 1, 



follows: 
where 



H Si = detBi/detB, 

T/(7l) T/(7fc) 

v si ■ ■ ■ y si 



n reads as 
(54) 



B 



V t 



(71) 



... Vs 



ilk) 
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is a k x k matrix, and Bi are obtained from B by replacing V s by HL — 6L- for all j — 1, ... , k; 



k 

H^H^ + J2H S V^\ i = l,...,r, i^ Sj for j = l,...,k, (55) 
i=i 

where H Si are given by (5$) . 

Proof. First of all, note that the above formulas for if, indeed constitute the Stackel transform, as 
Eq. (l54p is readily seen to imply the relations of the type d2J), namely 

k 

H^+J2H s y s ^ = a u i = l,...,k, (56) 

i=i 

cf. the discussion after (jSj). 

Now we only have to prove that the Hamiltonians Hi defined by (|54|) and ([515]) have the separation 
curve ( 153]) . As we have already mentioned above, the Stackel-equivalent n-tuples of separable commuting 
Hamiltonians share the separation relations provided (128]) holds. Therefore, in order to prove our claim it 
suffices to show that the separation curves (152]) and (153]) can be identified by virtue of (156]) . 

Indeed, upon plugging into (|52l) the relations 

n 

A-* = _ J> p (7,) A n- P) j = i f ... f k, (57) 

p=l 

that follow from (|50|) . collecting the coefficients at the powers of A, and taking into account (I5TI) . the 
separation curve (1521) can be rewritten as 

n 

"£Hf\ n -i=^(\,ii). (58) 

3=1 



On the other hand, plugging fl57|) into f[53|) and proceeding in a similar fashion as above, we obtain 

n / k \ 

~ E E^ (7j) + ^ • ■■ R n-s k (Fo) = V(A,A*)- (59) 

P =i \j=i / 



By virtue of relations (|56l) . which can be further rewritten as 

H$> = -j2n. i vt 1i) + a i , i = i 



j . . . j 



along with fl55l) . we find that the curves fl59l) and fl58l) are indeed identical, and hence so are the curves ([53 
and (1521). □ 



Remark 1 In fact the above argument can be inverted, that is, we can obtain the relations |37Jj (and hence 
(54\ ) ) and 153]) by requiring the curves / f5^j) and < f33|) to coincide and comparing the coefficients at the powers 
of X on the left-hand sides of these curves, or equivalently (by virtue of / f50j) ). of / f5Pj) and 155]) . 
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Proposition 7 The inverse of the k-parameter generalized Stackel transform (54\ ), 153]) has the form: 

H s .=detBi/detB, (60) 

where 

B -- 



T/ (n- Sl ) fr{n-s k ) 
V si ' ' ' V si 



v f 



(n-si) 



... V, 



(n-s k ) 



is a k x k matrix, and Bi are obtained from B by replacing Vs™ ^ by H^' — ctj for all j = 1, . . . , k; 



r(0) 



(n-sj) 



i = l,...,r, i^Sj for j = l,...,k, 



(61) 



where H Si are given by / TfiOj) and t^- are deformed separable potentials defined for all integer m except 
m = ji, i = 1, . . . , k, by means of the relations 



x m + J2 Vj^X^ + V p {m) \ n ~ p = 

j=l p=l,p^si,...,s fc 



(62) 



that must hold for A = Aj ; i — 1, . . . , n. 



The proof of this result is readily obtained from that of Proposition [6] using the fact that the inverse of 
the n-parameter generalized Stackel transform (|54j) . (|55l) is nothing but the dual of the latter (see Section [2] 
for the definition of duality). 

As a final remark, note that upon setting the parameters ctj and <Sj to zero for all % = 1, . . . , k the formulas 
f)54p and (1551) indeed relate the Hamiltonians Hi with the separation curve (1491) and the Hamiltonians Hi with 
the separation curve (jUJ). In this case we essentially recover the formulas from [7] relating the Hamiltonians 
from the seed class and from the so-called fc-hole deformation thereof (in our language, the deformed systems 
are precisely those having the separation curve PHI) ) up to a suitable renumeration of the Hamiltonians Hi. 



7 Examples 

As a simple illustration of the above results, consider the Hamiltonian systems on a four-dimensional phase 
space M = M. 4 with the coordinates G°i,P2, <?i, 02) anc ^ canonical Poisson structure. 

For our first example let k = 1, r = 2, s\ = 2, a% = a and ot\ = a. Consider the Hamiltonian 

rr 1 2 . 1 2 , ~ <£) 2 2 

#1 = oPi + 0P2 + P2 - 2a g l5 

which is Liouville integrable because it Poisson commutes with 

„ giP2 - O2P1 - 2agig 2 
"2 = • 

P2 

The above pair of commuting Hamiltonians was found by analogy with one of the models from [T9] . 
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The relation ([3]) in this case takes the form 

= a, 

P2 

whence 

5- ?iP2 - Q2P1 ~ ot-V-2 
2gig 2 

and therefore by virtue of (HI) we have 

~ q 2 1 +q 2 2 ~2aq 1 a(q 2 - aq x + g 2 ) 2 
1 = 9^7, PlP2 + P2- 

By Proposition [H ii) the relation {Hi,H 2 } = implies {Hi,H 2 } = 0, so Hi is Liouville integrable just 
like Hi. Interestingly enough, in this example the generalized Stackel transform sends the Hamiltonian Hi 
into a natural geodesic Hamiltonian Hi, but the metric associated with Hi is not flat and, moreover, has 
nonconstant scalar curvature unlike the metric associated with H\. 
By Proposition [3] the reciprocal transformation 

Z + a* ( o 1 (9i -2agi + gf> 2 \ 2gig 2 
V 92 / P2 

takes the equations of motion for Hi and if 2 , with the respective evolution parameters ti and t 2 , restricted 
onto the common level surface N& = {x G IR 4 |i?i(x, ai, a 2 ) = Si, H 2 (x, a\, a 2 ) = <S 2 } into the equations 
of motion for Hi and H 2 , with the respective evolution parameters t\ and t 2 , restricted onto the common 
level surface N a = {16 M 4 |i/i(x, ai, a 2 ) = «i, Hz(x, a 1( a 2 ) = a 2 }. It is easily seen that N a and N & indeed 
represent the same submanifold of M 4 . 

For the second example we set k = r = 2 and consider the (extended) Henon-Heiles system with the 
Hamiltonian 

TT 1 2 . 1 2 f 3, 9l9f\, 

#1 = gPi + 2 P2 ~ ttl \ qi + ~T~ J ~ a2qi1 

which Poisson commutes with 

2 / 92 , 9i 2 9lA _ 9! 



#2 = -92P1P2 - -91P2 - «i + —) ~ "2 j 
The separation curve for the system in question belongs to the seed class and reads 

ai\ 4 + a 2 X 2 + Hi\ + H 2 = \fi 2 /2. (63) 
The separation coordinates (A,, /ij), i = 1,2, are related to p's and g's by the formulas 



q 1 = \ 1 + A 2 , q 2 = 2V-AiA 



2- 



Pi 



Ai/ii A 2 /i 2 / — - — — / fii \l 2 
+ 7 P2 = v A^A 2 I — — + 



Ai — A 2 A 2 — Ai \Ai — A 2 A 2 — Ai 

Let si = 1, s 2 = 2, k = r = 2. Then fl5Tj) and fl54|) yield the following deformation of ifi and H 2 : 



~ 2 2 8 2(g|+4gj) 2 4 16 „ 

#1 = 2P1 ~ — ^1^2 4 P 2 2 a l + —"2, 

9i 92 92 9i 92 9i 92 92 



49i 2 + qj 2 4(g 2 2 + 2q 2 ) lQq\ + 12g 2 g 2 + q\ 2 

—r, 2~Pl 3 PlP2 + 4 P 2 

29i92 92 9i92 



, (92 2 + H) 8(g 2 + 2q 2 ) 

H 5 a i i "2- 

9i92 95 
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The corresponding separation curve reads (see Proposition [6]) 



if x A 4 + H 2 X 2 + 6nX + a 2 = A/i 2 /2. (64) 

Using Proposition [3] and proceeding in analogy with the previous example we readily find that the recip- 
rocal transformation fl26l) for the equations of motion restricted onto the appropriate Lagrangian manifolds 
in our case takes the form 

dh = U + d tl +(^ + ^f) dh, dt 2 = qldt, + q ldt 2 . 
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